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ABSTRACT 


The process of cold flat plane-strain rolling of 
rigid-plastic strain hardening materials is analyzed by Galerkln 
finite element technique coupled with the direct penalty method at 
steady state. Analysis is carried out both without and with strain 
hardening and the analytical results for normal pressure 
distribution, roll force and roll torque with strain hardening are 
compared with the experimental results. A detailed parametric 
study is carried out to show the effect of the process variables 
viz. reduction ratio, R/h^ and the coefficient of friction on the 
important rolling design parameters of normal pressure 
distribution, roll separating force and roll torque. The 
percentage under estimation of these rolling parameters with the 
assumption of perfectly plastic material is highlighted. A study 
of contours of normal strain -rate, equivalent strain., rate, 
deviatoric stress and equivalent strain for four selected sets of 
process variables is presented. It is concluded that the present 
model is capable of estimating the important rolling parameters 
for a range of process variables and easily accommodates 
experimentally measured variable coefficient of friction for- 
estimating more accurate results. 
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CHAPTER 1 


INTRODUCTION 

1.1. General Introduction 

Rolling is one of the most important industrial metalworking 
processes because of its extensive use in the production of a 
large variety of industrial steel products. Moreover, as is 
evident from the extent of published literature since the 
beginning of this century, rolling is apparently one of the most 
investigated of all the known metal forming processes. 

1.2. Mechanics of Rolling Process 

In reality as the material approaches the rolls it first 

undergoes elastic deformation till it reaches the yield point and 

a large plastic deformation zone follows. As the material exits 

the deformation zone there will bean elastic spring back. However, 

the magnitude of the elastic strains Is negligibly small compared 

to the large amounts of plastic deformation that the material 

undergoes during the rolling process. Typical orders of elastic 

-4 

and plastic strains usually found are 10 and 0.25 respectively. 
Therefore ignoring the elastic strains the material may be assumed 


to be 

rigid 

till it 

comes first in 

contact 

with the 

rolls 

and 

after 

it leaves the 

contact of the 

rolls 

Further 

during 

the 

deformation 

process 

in the deformation 

zone the 

strains 

are 


totally plastic. 

Many experimental and theoretical investigations have shown 
that the velocities, strains and stresses are not uniformly 
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distributed Across the sections. The velocity across the sections 
of the work in the unstrained zone remains uniform while in the 
deformation zone the outer layers of the metal towards the rolls 
move faster than the inner layers. Near the entry the outer layers 
have a less velocity than the roll periphery. But the velocity of 
the outer layers (represented by V ) of the work piece Increases 

fi 

as it moves towards the exit and eventually it becomes more than 

the roll peripheral velocity(represented by V ) as it reaches the 

exit. Although the transition occurs over a zone, because of the 

lack of exact data many researchers considered this as a point. 

The zone in which V is less than V is called forward slip region 

r s 

and in this region the friction force tends to retard the motion 
of the work piece. The point at which V^ls equal to is called 

the neutral point. 

Analytical expressions for the determination of the neutral 
point have been developed by some researchers. Notable among them 
is the Avltzur’s formula [ 1 ]: 


ot 


n 



tan 



1 


- 1 )- 



h 

< ln <h 


1 


) 



- 1 ) 


+ ----- - 1 H 

(2/y3)o- o 

The strain rates developed in the rolling process where the 
roll velocities are not very high are appreciably small. The 
entire work done by the tool in the plastic deformation of the 
workpiece is converted into heat leading to an increase in the 
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temperature of the workpiece. Strain-rates, temperatures and the 
strains that the material undergoes can affect its yield strength. 
However, the effect of the former two viz. strain-rate and 
temperature effects can be ignored in a cold rolling process since 
they are not substantially high whereas the strain hardening 
effect is substantial to be considered. 

Unlike most of the other metal forming processes where the 
tool remains stationary, in the flat rolling process the rolls 
will also be in movement thus the mathematical modelling of the 
process is difficult. The friction coefficient constantly changes 
from point to point along the arc of contact becoming zero at the 
neutral point. Many analytical expressions have been proposed in 
the literature to model the friction behaviour at the Interface in 
the rolling process but none has predicted its variation 
correctly. 

For a successful design of rolling equipment and the process 
proper the Important parameters to be considered are the normal 
pressure distribution along the interface, roll separating force 
and the roll torque. Analytical estimation of these parameters 
entails the analysis of the strainratea, strains in the workpiece 
during the deformation process. 

1.3. Previous Uork 

Many researchers carried out a variety of analytical and 
experimental investigations on the complex problem of rolling 
process . 

The pioneering work on the theory of rolling was done by Von 
Karman [ 2 ] in 1925, and then by Orowon [ 3] in 1943. As per the 
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review given In [ 5 ], von Kerman Introduced for the first time 
the slab method of analysis where the assumption of planes remain 
was made, while the later work came out to be one of the most 
comprehensive of all analytical studies where Orowon accounted for 
all the various phenomena of rolling process and developed a 
differential equation. Orowon could not give an analytical 
solution to it but he presented a graphical solution. Considering 
the complexity of Orowon’s equat ions , the later researchers, 
especially Bland & Ford [ 4 ] introduced some simplifying 

assumptions and developed analytical expressions. Alexandar[ 5 ] 

presented a comprehensive solution to Von Karman’s original 
differential equation without taking Into consideration the 
effects of inhomogeneity of deformation. Lahoti et al [ 6 ] 

combined the best of the theories then existing and solved for 
temperatures , roll separating force and roll torque with the aid of 
a computer-aided simulation and compared his results with the 
experiments in the literature. Numerous slip-line field and 
upper-bound solutions have also been proposed for the plain-strain 
rolling case. A few popular ones are the Alexander's [ 7 ] 

slip-line field solution for hot rolling, the slip-line field 
solution of Firbank etal [ 8 ]for cold rolling and the only known 

upper bound solution of Jhonson etal [ 9 3- The slip-line and 

upper bound techniques can give good solutions only to a limited 
number of cases and the material is invariably assumed to be 
perfectly plastic. The solution is also not unique in these cases. 

The advent of the application of the powerful technique of 
finite element method to metal forming processes is of recent ■, i 
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past and It Indeed paved a way to a great multitude of analytical 
aolutlona which till date have taken different turns towards the 
advancement in unraveling the problems involved in the rolling 
process, especially the neutral point phenomenon. The initial uses 
of FEH in metal forming problems were made by Zieklewicz etal 
[ 10 ] and Dawson[ 11 ] who treated the rolling problem in terms 

of vlscoplast iclty where the metal is considered to flow as an 
incompressible non-Newtonian fluid with prescribed boundary 
velocities. Dawson also Included the thermal analysis . G. J . Li . et 
al[ 12 ] have analyzed the problem of cold flat rolling by 

rigid-plastic finite element analysis by considering the friction 
relation developed by Chen and Kobayashl[ 13 ] . At the same time 
Mori etal [14) analyzed in the same direction but incorporated 
slight compressibility of the work material. Yhu et al [ 15 ] 

accounted for the variable coefficient of friction by using the 
experimental data in the literature and obtained a better 
correlation with the experimental results. Rlchelson [ 16 ] 

studied the effect of various friction models on the distribution 
of the normal and tangential tractions in cold plane-strain 
rolling using rigid-plastic FEM. However, the final important 
parameters of roll force and roll torque were not estimated in 
this work rendering it to be incomplete. An interesting case of 
asymmetric cold rolling was dealt with by Pletrzyk [ 17 ] where 

the problem of plane-strain rolling was analyzed with different 
angular velocities and different diameters of the rolls 
considered. The friction coefficient was used differently for 

the two rolls in the pair. Mori etal [18] tried to incorporate 
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the elastic deformation of the rolls by boundary element method 
while the strip was analyzed by rigid-plastic finite element 
method. It was Illustrated that the normal pressure Is lower with 
deformed rolls, than that with rigid rolls , in contrast to what 

was believed by many researchers. Some 3-D analyses have been 

conducted by researchers notable among them are [ 19-24 ] . 

A good number of experimental Investigations have also been 
conducted In rolling. Al-Salehl et al [ 25 ] has presented a 

comprehensive set of experimental results by measuring the roll 

pressure distribution in the roll gap during cold metal rolling 

under plane-strain conditions for aluminium, copper and mild 
steel. The difference in pressure distribution with different R/hl 
values was highlighted. It was shown that neutral point occurs in 
all cases of rolling very close to the exit. Roll force and roll 
torque values were measured and presented for aluminium and copper 
at various reduction ratios. Shlda etal [26] have measured the 
roll force and roll torque In case of cold rolling of mild steel. 
1.4. Objectives of the present thesis 

From the above study of the finite element literature it is 
clear that while the initial works ignored the strain-hardening 
effect, the later investigations except for [ 15 ], considered a 
frictional relation in which the so called friction factor 'm' is 
used which does not seem to have any rational basis. Also few 
works have estimated the error which results from the assumption 
of perfectly plastic material. Some researchers did not determine 
the practically useful quantities like the roll force and roll 
torque while many others did not compare their results with any 
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experimental work. 

In the present work the important design parameters o£ the 
rolling process viz. roll force, roll torque and normal pressure^ 
shear stress distributions are estimated for various cases 
involving three metals under the plane-strain assumption and by 
considering the effect of strain hardening. For treating friction 
at the work piece-roll interface and to model the neutral point 
phenomenon, the traditional Coulomb’s friction relation is used in 
which the average coefficient friction is taken from the same 
experimental works in the literature from which the analyzed cases 
are considered. The analytical results for roll force, roll torque 
and normal pressure distribution are compared with the 

experimental results. A detailed parametric study is carried out 
and the percentage error that is likely to accumulate with the 
ignorance of strain hardening effect is highlighted. 

Determination of the above mentioned design parameters of the 
rolling process needs the knowledge of stress and strain 
distributions in the work piece. These are estimated by treating 
the rolling problem as a plane-strain steady state problem using 
Eulerlan formulation with strain rate as the measure of 
deformation. The material is assumed to be rigid-plastic strain 
hardening. Temperature and strain-rate effects are ignored as 
these are usually small. For simplicity, the rolls are considered 
to be rigid. The domain includes sufficient portions of the inlet 
and exit zones to take the advantage of uniform conditions 
prevailing there. Galerkin weighed residual technique coupled with 
the direct penalty method is used for weak formulation. The domain 
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Is discretized Into small elements and biquadratic approximation 
for velocity components is used to convert the continuity and 
momentum equations into a set of nonlinear algebraic equations. 
These equations are solved by iteration. First a converged solution 
is obtained without the strain hardening effect. This solution is 
then used as the initial approximation to solve the system of 
equations for strain hardening problem. The knowledge of the 
equivalent strain required for updating the material properties at 
each Iteration is obtained by integrating the equivalent strain 
rate. The converged solution for nodal velocities is used to 
estimate secondary quantities viz. strain rat es , stresses , rol 1 
separating force, roll torque and normal pressure and the shear 
stress distributions. 

Our modelling is such that the roll velocity is not required 
in our non-dimensional formulation. 

1.5. Scope of the Present Work 

In this analysis, the material is assumed to be rigid plastic 
strain hardening which yields according to the von Mlses 
criterion. Although this assumption is valid in the deformation 
zone the elastic effects are significant at the entry and exit and 
therefore, should be included in the analysis. However, for 
simplicity, these effects are ignored in present analysis. The 
rolls are assumed to be rigid. Roll flattening has a propensity to 
decrease the roll forces since the effective reduction with 
elastically deformed rolls is lesser [ 18 ]• 

The effect of temperature and strain rate effects on the 
yield strength of the material are Ignored in this work because of 
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the two reasons, besides the fact that their inclusion renders the 
analysis to be too complex. The first reason is that it is 
strongly expected that whatever the direct effect of strain ™rate 
on the yield strength is present it is marginalized by the 
temperature softening. Secondly the roll speeds are not so high as 
to increase the strain rates to magnitudes. It Is observable that 
no previous work has considered both strain and strain rate 
effects simultaneously. 

The most important of all assumptions is that of average 
coefficient of friction. But in practice this is not the case 
where the coefficient of friction continuously varies along the 
Interface, becoming zero at the neutral point. 

The roll velocity can enter the problem only through the 
dependence of the yield stress on the strain rate or through the 
roll velocity dependent friction relation. But since we are not 
incorporating any of these two effects, there is no dependence of 
various parameters on the roll velocity in this formulation. 

The present analysis investigates into the applicability of 
the traditional Coloumb’s friction relation for the treatment of 
the friction at the workpiece-roll Interface in the problem of 
cold flat plane-strain rolling of plates where the width of the 
plate is considerably large compared to its thickness, at steady 
state. In this relation the frictional stress is a constant 
multiple of the normal pressure, this constant being the average 
coefficient of friction and in this analysis this has been 
obtained from Al-Salehi's[ 25 ] experimental results for aluminium 
and copper and from Shida etal [263 for steel. The deformation 
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of the plate is analyzed by the rigid-plastic finite element 
technique coupled with the direct penalty method and Eulerlan 
formulation. An endeavour la made in the analyala to llluatrate 
the effect of atraln hardening on the important parameter a of roll 
force, roll torque and the normal pressure dlatribution which are 
of utmoat importance In the mill design. The reaaonably good 
correlation between experimental reaulta [ 25,26 ] and theoretical 
reaulte over a range of input parametera ahowa that thia model haa 
certain validity and if extended with variable coefficient 
friction, non-rigid rolls and elastic entry and exit it can 
produce more accurate results. 

1.6. Plan of the Thesis 

The thesis is organized as follows. 

In the second chapter the mathematical modelling and the 
Galerkln finite element formulation of the problem of plane-strain 
cold flat rolling are presented. The application of the boundary 
equations and the solution procedure are also discussed in this 
chapter . 

In chapter 3 the results obtained by running the model for a 
range of cases are presented along with the comparison with the 
experimental results and a discussion on parametric study. A 
conclusion is drawn and suggestions for further work in the field 
are given. 
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CHAPTER 2 


MATHEMATICAL MODELLING OF COLD PLANE-STRAIN ROLLING PROCESS 

2.1. Introduction 

In this chapter the mathematical model for flat cold rolling of 
plates is developed. Considering that the width of the plate is 
more than ten times the initial thickness , the process is modelled 
as a steady state plane-strain problem. The constitutive equations 
for the mechanical behaviour of the work material are stated. The 
interaction between the work and its surroundings are represented 
by appropriate boundary conditions. Non-dimenslonalizatlon of 
governing relations and Galerkin formulation are presented at the 
end of this chapter. 

2.2. Governing equations 

The mechanical behavior of the metal is governed by the 
following two equations. 

(1) Conservation of mass or continuity equation : 

p + p div.V =0 (2.1) 

(2) Conservation of momentum or momentum equation : 

p V = (for zero body forces) (2.2) 

where the dot above V denotes the material time derivative. 
The constitutive relations , besides these governing equations , that 
are used to model the metal behavior are as follows. 

In flow formulation the measure of deformation is the rate 
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of strain tensor, which may mathematically be expressed as 


IT 1 j ] 

jt’V’V] °r 'ij= 1 Ciij* HP < 23 

In order to express the stress as a function of strain rate 
in a convenient form the stress and strain rate tensors can be 
divided into two parts_ 


g= -P 1+ | or = - pi5 + j (2.4) 

where p=-|tr a is the hydrostatic part, S is the deviatoric part 
and 6^ is the kronecker delta. Similarly, 


= i tr 


+ £ 


or 


ij 


^kk^ij 


+ £ 


ij 


(2.5) 


where £ is the deviatoric part of the strain — rate tensor. In 
plastic deformation, since there is no change in volume the 
hydrostatic part 1s not related to the deformation. Another 
consequence of volume constancy is that the hydrostatic part of 
the strain rate is zero and its deviatoric part is the strain rate 
tensor itself. Then the deviatoric parts of stress and strain rate 
tensors can be related as : 


5 = 2/j£ or S i j = 2 jj £ ^ ^ (2.6) 

where 2(j is the proportionality constant. Further if we define 
the second invariant of deviatoric stress as 


§ "^i s ij s ij 


and second invariant of strain rate 


(also known as 


(2.7) 

equivalent 
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strain rate ) as 



2 *ij*ij 


then it follows that 


( 2 . 8 ) 


S=yj£. (2.9) 

For metals which yield according to von Mlses criterion. 


a 

S= - y 


( 2 . 10 ) 


and yj for such metals comes out to be 

w- "y /y-? I (2-11) 

where o- is the yield strength of the metal. In general & may be 

y y 

a function of equivalent strain , temperature and equivalent strain 
rate. However, if the deformation process is such that the 
temperature and visco-plastic effects on the mechanical properties 
can be neglected then it is a function of only equivalent strain. 

For a strain hardening metal the yield stress is related to 
equivalent strain by an equation. 


=©■ (1 + K£) n (2.12) 

y o v 

where <r is yield stress of the metal at zero plastic strain and K 
and n are metal dependent coefficients determined from 
experiments. The equivalent strain is obtained by time integration 
of equivalent strain rate as given below! 

t 

£=J J dt (2.13) 

0 

2.3. Equations for steady-state two-dimensional problem 
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The work 

is assumed to 

be 

very wide 

as compared 

to 

its 

thickness and 

consequently 

the 

effect of 

the spread 

may 

be 


neglected and the problem may be considered as plane-strain. The 
rolls are assumed to be rigid and the metal is assumed to be 
rigid-plastic strain hardening yielding according to von Mlses 
criterion. The temperature and visco-plastic effects are neglected 
because it has been observed that the strain rates in the cases of 
cold rolling process considered are ignorably small. Only half of 
the domain is considered for analysis because of symmetry, see 

figure 2.1, and the rectangular cartesian coordinate system as 

I 

depicted in figure 2.2 is used. 

For plane-strain problem the velocity vector has only two 
components , 

V=ui+v j (2.14) 

while the stress and strain-rate tensors have the following form: 



<y 

XX 

<y 

xy 

£■ = 

jt 

XX 

xy 


<y 

L y* 

<y 

yy 

— - 

£ 

yx 

£ 

yy J 


c 


The process is analyzed at steady state and therefore the 
governing equations as reduced to 2-D are: 


<9u 

Sx 



(2.16) 


c du^ an 

C v iy 3 


a as as 

+ XX + xy 

ax ax ay 


r av , av •. 

C u i;* v iy ) 


. . as as 

#2 + __yx + yy 

a y ax ay 


(2.17) 


(2.18) 


where 
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S =2fJ£ =2p" 

XX xx Sx 


S =2>J£ =2m?- 

yy yy *y 


S = 2fj£ 
xy xy 


=fJ ( 5y + ii ) 


(2.19) 



2.4. Boundary conditions 

The boundary conditions specified on any boundary can be 
classified as essential boundary conditions i.e., the 
specification of the velocity vector V, natural boundary 
conditions i.e., the specification of the traction vector t=^.n or 
mixed boundary conditions i . e ., sped f icat ion of one velocity 
component and one traction component. For the problem under 
consideration the boundary conditions are as £ ollows ( f igur e 2.3). 

(1) Entry and exit boundaries (AF&DE): 

The control volume is so selected that its entry and exit 
boundaries are sufficiently far away from the either side of the 
deformation zone. This ensures a uniform variation of the velocity 
before and after deformation and allows the following 
boundarycondit ions for consideration on these regions! 


li 

3 

and 

o 

11 

> 

on side 

DE 

u=U 1 

and 

o 

II 

> 

, on side 

AF 


(2.23) 


i.e., all the material particles on these boundaries have the 
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prescribed velocities and at entry and exit respectively, 
and the velocities of these particles in Y-dlrectlon are zero. 

(2) The top free surfaces ( AB&CD) : 

Since these surfaces are traction free, 

t =0 and t =0 
x y 

Also since the flow is only in the X-dir ect ion , v=0 . But only 
one condition out of the two (t =0,v=0) has to be chosen. So the 

y 

final choice of the boundary conditions is 

t =0 and v=0 (2.24) 

(3) The axis of symmetry(FE) : 

Because of symmetry, v and t are zero on this region of the 
boundary. Thus 

v=0 and t x =0 on the axis of symmetry. (2.25) 

(4) Uork-roll interface (B€): 

The component of velocity in the direction normal to the 

roll-plate interface at any point on the Interface is zero. 

u =0 or u+v tanot=0 (2.26) 

n 

where u and v are the velocities in the directions x and y at any 
point on the interface and a is the angle made by the unit normal 
to the interface at that point with the y-axis . At the first 
contact point on the interface towards entry the angle ot is equal 
to the angle of bite. 

The complex and important part of the analysis of rolling 
process is the treatment of the work-roll interface friction and 
modelling of the neutral point phenomenon. Many researchers have 
tried to tackle this problem by using rather arbitrary functional 
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relations for shear stress in such a way that the position of 
neutral point is constrained to be within the inlet and exit of 
the interface, in accordance with the experimental findings. 
According to Al-Salehl etal,[ 25 ] the neutral point, such that 
the coefficient of friction falls to zero and a shear stress 
reversal occurs, is present in every case and close to the exit 
plane. However, none of the works ha^ tried to incorporate the 
experimentally measured coefficient of friction either directly or 
indirectly. Such an endeavor is made in this analysis and thus the 
following relation is used: 

I n l < 2 - 27 > 

where f is the coefficient of friction. In this analysis an 

♦ * 

average value of fls used assuming that it remains constant 
throughout the deformation process. These values for individual 
cases considered are taken as reported by Al-Salehi etal . 

Further the shear stress in the eq . (2.27) is subject to the 
constraint 

o a 

if t > - y then t = - y 

S V3 S V3 

The method of application of this boundary condition is explained 
later in the next chapter in detail. 

2.5. Non-dimensionalization 

Non-dimensionalization is usually done to avoid 
ill-conditioning and numerical difficulties. The 

non-dimensionalization of all physical quantities using 
characteristic dimensional quantities is presented in the 
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following lines: 


rw 

x = 


X 




U = “ 


u 


U 


2 


v*-g- (2.28) 

2 


I f we use 


_o 

^o' y3C2U~/h") ’ 


(2.29) 


to non-dimenslonalize fj then 


the expression for fj may- 


be 


written 


as , 



(2 . 30) 


where the non-dimensionalized strain rate invariant is given by 



(2.31) 


For a perfectly plastic material 


P = V* 


(2.32) 


The non-dlmensionalizatlon of the other quantities is as 
follows : 


P= VV S 2 = "%/ 95 

S S 

~ xy xy 

S xy = ^;2U 2 7h2 = _ ^;7V3“” 


xx 


xx 


“xx ^ c 2U 2 /h 2 * 0 / y3 

c q 

s =___ _yy— _YY— 
yy ^ 0 2 V h 2 °'o/ y3 


(2.33) 


The non-dimensionalization of the continuity and momentum 
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equations ia obtained by substituting the eqs.(2.28) to (2.33) in 
eqs .(2.16) to (2.18) . 


-z + ~b° 

Sx Sy 


(2.34) 


pu h . ~ ~ % ~ , SS SS . 

-H 1+ -£ _f __** + __*Y 1 =0 

Sx Sy ' Sx ^ Sx Sy * 


similarly, 
^U 0 h 


f~ Sv ~ 
|U --+V 

Sv 

V 4i -f 

. SS SS . 

+ __YY 1 

^ Sx 

Sy 

j ay ' 

■ Sx Sy ' 


)=° 


(2.35) 


(2.36) 


Since the Reynolds number 



is found to be typically of 


“8 

the order of 10 for all the cases considered, the inertial terms 
may be ignored. Thus the summarised non-dimensional forms of the 
continuity and momentum equations are: 


IV «v 



Sx Sy 


(2.37) 


- + 


Sx 


, SS SS 
f __xx + xy 

^ Sx Sy 



(2 . 38) 


- *2 + f --VX+--YV 1=o (2.39) 

Sy ^ Sx Sy ' 


2.6. Galerkin Finite Element Formulation 

2.6.1. Gal erkin Ueak Formulation with Direct Penalty Method: 
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In this analysis the Galerkln weighed residual method 
coupled with direct penalization Is used for the finite element 
formulation of the governing equations. In the mixed 
or pressure-velocity formulation where pressure is also a nodal 
variable, the coefficient matrix becomes a non-posit ive-def lnlte 
matrix because of the zeroes appearing along the principal 
diagonal, therefore entailing a solution method employing complete 
pivoting. Also since the pressure does not occur as a degree of 
freedom at every node, the computer implementation of assembly 
becomes more complicated [ 27 ]. To overcome these problems and 
because of some of the other advantages with direct penalty 
formulation with reduced integration [ 28 ], this method is 
adopted in the present analysis. 

In direct penalty formulation the pressure In the momentum 
equations is eliminated using the following relationship 


where 



(2.40) 


and X, which Is called the penalty factor, is a very large 
number. In the limit as X -+ co , equation (2.40) reduces to equation 
(2. 37). The value of X is usually chosen in such a way that 


the upper limit being dependent on the computer. 

The weak form is derived In this sub-section and the finite 
element formulation is presented in the next sub-section. Let u 
and v be the functions which satisfy all the essential boundary 
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conditions exactly. Then u and v constitute a weak solution if the 
following integral equation is satisfied. 


I 


w + *i w 

Sx u Sy v 



as 

as 

„ , as 

as . ->1 


XX 

+ __*Y 

K - y - 

+ 1) 

dA =0 

ax 

Sy 

* ax 

Sy 



(2.41) 


where w^ and are weight functions which satisfy the homogeneous 
versions of the boundary conditions and 'A' represents the area of 
the domain. 

Performing the above integration by parts the following weak 
form can be obtained: 


I 


aw 

S\j 

, S\j 


r 4w av i 

~ ^ w . 

i 

! 

\C 

f- 

1 

1 

\< 

J s __H 

4- Q 

1 

1 

h 

i 

i 

\< 

4. Q 

ax 

Sy J 

1 XX ~ 

v ax 

xy 

~ t ~ 
[ay a x J 

* O 

yy *y 


dA 


I I fw S n + w S n +w S n +w S n 1- fw pn +w pn 1 dT =0 
J u xx x u xy y v yx x v yy yJ u x v^ yJ j 

T (2.42) 

where r represents the boundary of the entire domain under 

\ 

consideration . 

Substituting for p from eq . (2.40) and for form 

eq. (2.19) only in the area integral, the relation in eq . (2.4) 
only in the boundary integral and using the relation 


t y- ,L yx yy J t y J 


(2 . 43) 
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we obtain the final weak form as 


dv a\j 


\ Sf^lf 

J ^ ax ay J ^ ax av J L 


Sx a y 


a\j ~ 

| JC +jC 

XX xy 



- ft w +t u ] dr =0 
J L x u y vJ 


(2.44) 


2.6.2 Finite Element Approximation: 

In contrast to the traditional Gal erkin-Ueighed Residual 
Method where a global approximation is chosen for the field 
variable, in finite element method the domain is discretized into 
small elements and the above integral is evaluated over each 
element after suitably selecting the approximation for the field 
variable over each element. Then the elemental expressions are 
assembled to yield the global FEM expression. 

The eq. (2.44) suggests that the convergence criteria will be 

satisfied if u and v are chosen to be bilinear functions of x and 

*“ v 

y. So one can use bilinear interpolation for u and v. 
However , higher polynomials give better results. Ue choose 
nine-noded Lagrangian element and corresponding biquadratic nodal 
shape functions for u and v, figure 2.4* 

The approximation for u and v Is as follows: 
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(2.45) 



*••■.[ 5 ]- [ » v ][ n 

In Galerkln formulation the weight functions are approximated 
using the shape functions as that for velocity. 


[“]={3“} '[ N v]C ue ] C2.46) 

I 

The vectors [V 6 ] and [“*] respectively contain the velocities 
at the nodes of elements and their corresponding weight functions. 

The expressions for the shape functions for nine noded 
Lagrangian quadratic rectangular element are given by: 


N i”i P 2 -0 C« 2 -") > 
Vi P 2 -OC« 2 *«) = 

»7*2 C 1 -' 2 )^ 2 *") ' 


Vi (« 2+ «)0» 2 -»D 

Vl (l-? 2 )(" 2 -«) 

vl p 2 -«)ci-» 2 d 


vl 2+ <) C” 2 *”] 


vl c« 2+ o c 1- ” 2 ) h 2 - 47) 


V C 1 "* 2 ) ( 1_T > 2 ) 


where £ and r> are the local coordinates and are related to the 
global coordinates x and y with subparametr ic formulation as 
follows : 
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(2 . 48) 


4 

x=£ x jL N i 


where 




N 3 = i( 1_<: ) ( 1+n ) 


4 


y=E v i N i 

L*=i 


ViC^OC 1 -*) 

(2.49) 

ViO'OC 1 *") 


The boundary approximation, figure 2.5, is made consistent 

with the area approximation as follows. 



[ "J [ <£*] 


(2.50) 


where ["J is the vector of Lagrangian quadratic shape functions 
and is the vector of local boundary weight functions for a 

typical - boundary element. Further we approximate t x and t 

as , 



(2.51) 
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t> & 

where -{t^ } is the vector of local boundary tractions for a 
typical boundary element. 

2.7. Finite Element Equations in Local Variables 

The following definitions of some vectors should precede the 
presentation of elemental expressions and their assembly, for 
compactness and ease of understanding. 

The strain-rate vector is defines as follows. 


or further 



(2.52) 


(2.53) 
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and 


where 


aw 

u 

~Zi 

a x 

S\j 


ay 

. r dv aw . 

if— “+_ y 

V2 U 




* [ ® ] E U6 1 


(2.54) 


[ » ] 


1 

V2 


SN^ 

IV 

ax 

0 

Wj 

IV 

ay 


o 

aj^ 

ay 

1 ,N 1 
y2 


aN, 

i 

Sx 

0 


aN 

1 2 

72 ~ 

ay 


0 

f_2 

ay 

1 ™2 


V2 


ax 


aN { 

ax 

o 


aN 

1 __9 

y5 ay 


0 

aN ? 

ay 
i aN 5 

^ ax 


(2.55) 

from 


Substituting the the area elemental approximation 
eq . (2.45) to eq. (2.49) in the area terms of eq. (2.44) and 

noting that 


au + av 
ax ay 


= [»] T [B][V e ] 


and 


a w aw 

-7* “ T -[ »•][*][»] 

ax ay 


where 


[ "] 


1 

1 

0 


.we obtain the elemental area term in matrix form as 
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J [£ {[ Ot B] T [»][» ] T [ B][? e ]V 

A e 

{z2[ U*] T [ B ] T t B ][ V*]jJdA e (2.56) 

Further substituting the boundary approximation from eqs . (2.50) 

and (2.51) in the boundary term of eq . (2.44) we obtain the 

elemental boundary term as 


J[“b e ] T [»J I [N b 3< ^ <“- e 

,-,e 


(2.57) 


To convert the above integral equations from the global 
coordinates to local coordinates ,x and y are to be changed to £ 
and n using the following relation. 

+ 1 +1 

J( ) dx dy = J | ( )|J| d? dr, (2.58) 

-1 -1 


where |J|,the jacobian,is the determinant of the jacobian matrix, 


dx dx 

d? dY? 

dy dy 

d? dr, 

Thus the eq. (2.56) and eq. (2.57) can be written in terms of 
natural coordinates as 
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+ 1 +1 


1 1 [M [ue][B]I[ ” 3[ " ]ItB]t ^} |j| « d " + 

-1 -1 


and 


+ 1 +1 

J | {tf[»*] I [B] T [B][S*]}]|J| d? d„ 

-1 'I 


(2.59) 


+ 1 

J[ u b e ] T [ N b ] T [ |, J<*b e H J bK f 2 - 60 

-1 

respectively where |J| is the jacobian for the area terms and JJ^j 
is the jacobian for the boundary term and £ is the natural 
coordinate for the boundary linear elements. 

2.8 Global Assembly 

The final finite element equation is obtained by 
assembling the elemental area and boundary matrices into 
global area and boundary terms , respectively , and then equating as 
follows . 
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nel em 

£ 

e=l 


+ 1 +1 

lj J{[ u6 ][ B ] T [ ro ][ » ] T [ B ][ de d 77 + 

-1 -1 

+ 1 +1 

J J { 2 ^[ “*] T [ B ] T [ b ] [ v*j}|j| dir 

-1 -1 

+ 1 

nbelemT r 


Further it can alao be written aa 


nel em 
£ 

e=l 


[“Jft **][*•] 


nbelem 

£ 

be=l 




where 


[**J- [*V*l > e ] 2 

+ 1 +1 

[ | J {*[ B ] T [ B ]}|J| d? d„ 

-1 -1 


+ 1 +1 

[**3 2 - J J{[ b ] T [”]["] T [b]||J| « *> 

-1 -1 

+ 1 

[ { l> J{[ N b3 T [ N bP ? ^}l J bl « 

-1 


b« dC j 
( 2 . 61 ) 


( 2 . 62 ) 


( 2 . 63 ) 


( 2 . 64 ) 


( 2 . 65 ) 


( 2 . 66 ) 
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The elemental area matrix £k®J is prepared by evaluating the 
integral by reduced integration technique where [k®] is 
integrated by 3*3 gauss quadrature and [k®] 2 by 2*2 gauss 
quadrature. This ensures that £k J ^ i® singular so that even for 
fairly large values of X,the penalty factor, a non-trivial solution 
may be obtained. 

The assembly of the elemental area matrices into the global 
area matrix is done by transferring the elements corresponding to 
a local degree of freedom in each elemental area matrix to 
positions of corresponding global degrees of freedom in the global 
area matrix. Similar procedure is followed for the assembly of 
global boundary matrix. The assembled finite element equation 
may be written as 

[U] T [K][V] = [U] T [F] (2.67) 

where [ U J contains global nodal weight functions, J" K ] is the 
global stiffness matrix, [ V ] contains the global nodal velocities 
and [ F ] is the global load matrix. 

Since [ U ] is arbitrary we can eliminate it from both sides 
and thus we obtain the final general equation of FEM, 


[ K ][ V ] = [ F ] (2.68) 

2.9. Application of Boundary conditions 

The application of natural and essential boundary conditions 
is done as follows. 

The natural boundary condition on the work-roll interface can 
be applied without requiring to evaluate the matrix [f£®] in the 
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following manner (figure 2.6). 

Let us consider a typical boundary element on the work 

piece-roll interface. The elemental boundary matrix for this 
element would be as given by eq. (2.66). At a middle node 'l' of 
this boundary element j, whose corresponding global node 

number is 'k^ we can write the contribution of this node to the 
local boundary matrix and there by to the global load matrix as 
follows . 

+ 1 

f ^ 

[ F Lic-rf f b Li-r J [ k J t [ N J) J b l 1 21 -i and 

>■ j 

-1 

+ 1 

C F ] 2 k “[ f b] 2 i “ ( |t "J T E » J I 21 

-1 

where [ F ^k-l anc * [ ^ 1 2k are e i eroents of the global load 

matrix containing the odd and even net force contributions from 
the global boundary node 'k’ t respectively. Now for this node the 
relation in equation eq. (2.27) can be written as 
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Converting the (t y ) k and (t ) fc Into matri etc analogous to, but 1 - 
boundary elemental traction matrix and multiplying both aides with 
+ 1 

T 

Jt N £][ N £l l J b f d < where [ wjj = [ Nj N 2 N 3 ] (2.68) 

-1 

and replacing (t' x ) fc and (t ) R with i *b> 2 i-l and { *&2i 
respectively, the following expression is obtained_ 

+ 1 

[j[ N b]L N J l J J d < < ? b> 1 2i-l [cosa k -fk s sin« k ] 

-1 

+ 1 

= [j[ N b3[ N J l J b l ^ < ^b>J 2 i [ einot k +£k s cosot k3 

-i 

This gives the contribution of one node of the element j'. 

Since ot does not vary much from one local node to 

other in an element , siml l&r evaluation for other two nodes and 
then summation over all three nodes of the element gives 

C £ b] 2i-l [ cos<2 ‘k" fk s sina k] = C £ i 2i[ 8inot k +£k s COSa k] • 

Further 

t F 3,2k- 1 [ COaC, k^ fk e ein “k] = t F 32k [ slna k* £1C » C °“'k3 

Thus instead of evaluating [ F ]2k-l and C ^ 3 2k £or any 
typical global interface node k, we perform the following the matrix 
operations in order to apply this natural boundary condition. 

For each boundary node 'k 1 on the interface, 

(1) replace each element of (2k-l)th row of global stiffness 
matrix by: J 
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[ element of (2k-l)th row • ( cosot, -fk sinot )- 

K s k y 

element of (2k)th row in the same column (sinot, +fk coaot. )1 

k s k J 

(ii) replace (2k-l)th element in global load matrix by 2 ero. 

The boundary condition expressed in eq. (2.26) at a typical 
global node V la applied through the following steps. 

(1) Replace all the elements of (2k)th row of the global 
stiffness matrix and that in the global load matrix by zeroes. 

(ii) Substitute ' 1 ' in the elemental position of (2k)th row and 
(2k-l)th column, and the value of tanot , evaluated with ot at node i^ 
in the elemental position of (2k)th row and (2k)th column of the 
global stiffness matrix. 

2.10. Formulation for Strain Hardening 

The equivalent strain may be obtained by the time integration 
of equivalent strain rates as, 


; ~vlr 


£ dt 


(2.72) 


where * and i> are the equivalent strain and strain rates 
respectively and the factor y| is introduced so that the 
equivalent strain is equivalent to the axial strain in a tensile 
test. Differentiation of the above equation leads to 


d£ It 
dt ~Y3 * 


(2.73) 


where ^-is the material time derivative as it is in Eulerian 
dt 

formulation. For 2-D steady state case, as ours, the above equation 
can be written in non-dimensionalized form as 


U-“ 


+ 



(2.74) 
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This equation has to be solved along with the condition ^=0 at the 
inlet . 

Gal erkin formulation of the above equation over an element 
gives the following integral: 


dx dy=o 
A 6 


(2.75) 


where w g la the weight function. 

The Integral may be expreaaed In matrix form aa given below. 


I [“< 


t u v] 


dx 

&£ 

*y 


- w 


a V3 


I *] 


dx dy=0 


(2.76) 


In the finite element approximation for equivalent atraina and 
the weight function the aame ahape functiona aa uaed in the 
velocity approximation are uaed. 


*=E N l = i n'> t < 
1 


V a=f "iW^ 


N l H 


a 


(2.77) 


where { N }- is the 1-D array of quadratic shape functions as 
defined below, 

i N , K = [N 1 N 2 N p ] (2.78) 

and the differentiation of atraina in equation (2.77) leada to 
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(2.79) 


SjC 

A/ 

aZ 

Sy 


[ « ]< 


where [ H ^ 


is defined &s given below. 




*3 


an 

— 

0 


0 ... 

— 

<9x 


ax 


ax 

0 


0 

f*2 

0 


ay 


«*y 



0 

*n 9 

Sy 

(2.80) 


Substituting the approximation for equivalent strain , weight 
function and velocities in eqs. (2.77) and (2.45) we obtain the 
following Integral equation. 


nelem 

r 

e=l 



<T[ b1 ] t [ v *] t [ » v ] t [ * b 


dx dy- 


yij ( [ v s e ] T [ I,1 ] T *] dx d y 


= 0 


(2.81) 


To facilitate the numerical integration, the variables are changed 
to the natural coordinates. The final expression can be obtained 
as 



nelem 


E 

e=l 


([ * > T = 


nelem 


A [' ’1) 


(2.82) 
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where 


+ 1 +1 

t ' | J {[ Nl ] T [ n T C » V ] T [ » ]} |J| « d„ (2.83) 

-1 'I 

la the elemental coefficient matrix for the evaluation of 
equivalent strains and 

+ 1 +1 

[ r e ] = y! J J|[ N 1 ] 1 ! j |J| d? dy> (2.84) 

-1 -1 

T 

is the elemental right side vector. Since £ U j is arbitrary, we 
obtain the following equation in terms of global matrices, 

[ K s ]< ^ ^ = C R ] (2.85) 

where the is the global coefficient matrix, •{£} is the global 

nodal equivalent strains matrix and [*] is the global right hand 
side, for the evaluation of equivalent strains. In equation (2.84) 

£ is the non-dimensional equivalent strain-rate which may directly 
be obtained at each gauss-point from the velocity field of 
previous iteration. 

2.11. Implementation 

Uhile correct modelling of the physical problem into a 
mathematical equation forms the core of the entire analysis, 
selection of degree of refinement of mesh, value of the penalty 
factor, the convergence criteria, the threshold value of ju and 
other mesh related parameters decides the achievement of 
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reasonable results. It has been experienced that the judgment 
about the later mentioned selection of various parameters needs a 
vide knowledge of capabilities and liabilities of computers as 
well it entails a good deal of numerical experiments before a 
particular set of values is used for final results. A brief 
discussion is presented in the following lines about some of these 
factors . 

The length of the inlet and exit zones is usually taken as 
follows : 


l^=nh^/2 and l^snh^/2 

where 1^ and 1 are the inlet and exit zone lengths respectively, 
'n' is a multiplying factor and are the inlet and exit 

thicknesses of the plate. The value for ‘n* has to be selected such 
that uniform conditions prevail at the beginning of the inlet zone 
and end of the exit, A value of 6 has been selected in this 
analysis after conducting a few numerical experiments. 

The value of penalty factor X is usually selected as per the 
directions given elsewhere in this chapter. For this analysis X is 
found as follows: 

X =X|U 

o 

~ 2 5 

where X is a numerical value in the range 10 to 10 . Although the 

pressure remains fairly unaltered for a reasonably wide range of 

values of X, very high values may render the solution to be 

spurious. A value of X=10 has been chosen in this analysis and 

used constantly for all the cases. 
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I** this analysis o£ rolling the neutral point Is not known a 
priori. Therefore, for the first iteration an approximate neutral 
point is determined by the following Avitzur’s analytical formula: 



& xb ~ °xf 1 

C2/V3)^; _ ] 


where 

ot = the neutral angle 
n 

R = roll radius 

h^ h 2 = the initial and final thicknesses of the plate 
m= the friction factor 

o'xb - the back and front tensions respectively 

a = the initial yield stress of the metal 
o 

Thus the relative signs k g at each interface node is 
determined as follows: 
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This approximation is used only for the first iteration and the 
second iteration onwards the relative signs are determined from 
the computations of t and t fi from the solution of 
the immediately preceding iteration (the computation of t 
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converged solution 


and 't la explained later). 

The strategy followed for obtaining the 
with and without strain hardening is as follows. First a converged 
solution (nodal velocities) without strain hardening has been 
obtained. For this portion of the analysis the p is simply as 
given in the eq. (2.32), This solution is then used as the initial 
solution for the further analysis and the yield strength is 
updated by evaluating strains from immediately preceding solution 
in every successive iteration. Thus the value of p is evaluated 
from the eq . (2.31). 

Preliminary runs have been conducted with coarser mesh using 
various solution algorithms viz. gauss elimination method, IMSL 
routine with and without iterative solving, NAG routine etc. on the 
HP super mini computer to check the validity of our gauss 
elimination subroutine. Then final runs with refined mesh are 
conducted on CONVEX mini super computer to harness the high speed 
parallel computing using this solution routine. Because of the 
assumption of rigid portions in the inlet and exit zones, the high 
value of p renders the coefficient matrix to be ill-conditioned, 
especially during strain hardening because the p is updated. To 
counter this problem, p is to be restricted to a threshold 
minimum, which is equal to the reciprocal of set off value of 
non-dimensional equivalent strain-rate. 

The convergence criteria used for velocity are as follows. 
All values of velocities less than 10~ 4 are considered to be equal 
to zero and are not checked for convergence. A percentage 

_ M _ 2 

difference leas than 10 for u-velocitles and 10 for 
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vvtlocitlsfl bfitwseft solutions of two successive Iterations is 
considered to mark the convergence of the solution , i.e., nodal 

velocities. This is done because the v-velocities are about two 
orders magnitude smaller than u-velocities . 

2.12. Evaluation of the Secondary Quantities 
As explained above the solution to the problem is obtained in the 
form of nodal velocities. Thereby the secondary quantities viz. 
strain rates, deviatoric stresses, hydrostatic pressure and 
interface boundary tractions are evaluated in each iteration to 
determine the non-linear parameters like £,k etc., as well as at 

£S 

the end from the final converged solution for parametric study. 
The details of these calculations are explained in the following 
lines briefly. 

(1) The Hydrostatic Pressure at Gauss Points: 

The pressure at the gauss points Is post-evaluated from the 
Lagranglan multiplier as follows: 

p=-\[ 1 = - X i m > T [ B ]-{ V% 

^ dx Sy J 

This is evaluated at 2x2 gauss points. 

(ii) Interface Boundary Tractions: 

Just like pressure the strain rates and the total stresses are 
also evaluated at 2x2 gauss points. The reciprocal of the 
non-dimensional strain-rate invariant helps in calculation of fj . 
The stresses at gauss points are linearly extrapolated on to the 
boundary nodes in a few steps. Then the tangential and normal 
tractions are evaluated at each interface node using the value of 
angle oi at that node. This is done in each iteration to find the 
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values of fj & as well as at the end for parametric study. 

(ill) The Roll Separating Force and Roll Torque: 

The roll separating force and torque are determined after 
obtaining a converged solution as follows. 

(a) Roll separating force: 

The roll separating force per unit width of the plate is 
evaluated as the Integral of the vertical traction component (t ) 

y 

over the arc of contact, expressed in terms of MN/m. Thus 

1 

s 

Roll separating force F =J t ds 

y y 
o 

where 1^ is the total length of the arc of contact. However, in 
this analysis, this integral is evaluated numerically using 2 
point gauss quadrature_ 

2 

F y:5 1 (< Vi l J l V 

where t is the normal traction at the gauss point *i,jJ| is the 
y 

Jacobian of the transformation from the boundary coordinates to 
the natural coordinates and w^ is the weight function. 

(b) Roll torque: 

The roll torque is evaluated by numerically integrating, by 2 
point gauss quadrature , the elemental torque along the arc of 
contact as given below: 
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Roll torque T r =J Rt fl ds =R £ (( t g ) x | J | w.] 

0 i = 1 

where " R " la the roll radius," t ” la the tangential traction 
at the gauss point and |J| and v i are as defined above. 

(iv) The Contours: 

Contours are plotted to s how the distribution of certain 

A/ 

secondary parameters across the domain, namely for 2 ,£ and S. The 
contours for £ are plotted using the data at the nodes and for the 
later two, the plotting is done using the values evaluated at 2*2 
gauss points of the elements from the converged solution. 
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V = o 



Fig. 2.3 The Boundary o£ the Domain and the Boundary Conditions 
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Fig. 2.4 A Typical Finite Element 
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Boundary Element tor Velocity 


Boundary Element tor 
Coordinate Tronstormation 


Fig. 2.5 A Typical Boundary Element 
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Fig. 2.6 The Boundary Tractions 




Fig. 2.7 The Finite Element Mesh 
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CHAPTER 3 


RESULTS AND DISCUSSION 

The finite element modelling of the physical problem of cold 
flat rolling developed in this work has been applied to a number 
of cases involving three metals and various non-dimensional input 
variables to illustrate its appl icabi 1 ity . The metals considered 
are aluminium , mi Id steel and copper whose stress-strain relations 
are shown in table 3.1. The materials , dimensions and the 

strain-strain relations are considered in such a way that they 
are consistent with the cases considered in the works 

[12,15,25,26 ] so that the comparison is facilitated. In this work 
only comparisons with the experimental results are presented and 
rest of the emphasis is put on presenting a detailed parametric 
study and the effect of strain hardening. Comparison with the 
other analytical results is not given. Thus, in this chapter, the 
analytically obtained results are presented in four sections. In 
the first section the analytical results are compared with the 
experimental results of Al-Salehi et al [ 25 ] in case of 

aluminium and copper and Shida et al [ 26 ] in case of steel 
to validate the modelling and computations. In the second 
section, the parametric study, the effect of the three important 
rolling process variables viz. reduction ratio, the homogeneity 
factor (R/hl) and the coefficient of friction on the rolling 
parameters viz. normal pressure and shear stress distributions 
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alone the arc of contact, roll separating force and roll torque 
are discussed. In the third section of this chapter, the error 
that can enter by ignoring the strain hardening effect is brought 
into light by comparing all the important design parameters 
without strain hardening and with strain hardening. The fourth 
section marks the end of the presentation of the results where a 
study of contours of normal strain —rates , equivalent strain 
rates , deviatoric stresses and equivalent strain is presented. In 
the penultimate section of this chapter conclusions are drawn as 
a post-analysis and evaluation of the work. In the last section 
suggestions for further research are included. 

3.1. Validation 

The validation of the results of the variation of the roll 
force (per unit width) with reduction ratio is presented in 
figures 3.1 to 3.6. The figures 3.7 to 3.12 show the comparison 
of the variation of the roll torque (per unit width) with 
reduction ratio. Finally the comparison of the distribution of 
pressure along the arc of contact is presented in figures 3.13 to 
3.16. 

It is clear from the figs.( 3.1 ) to ( 3.6 ) that there is a 
good correlation between experimental and theoretical results of 
the roll force in case of steel and aluminium and to some extent 
that for copper. However, at higher reductions the deviation 
between them increases except for the case of R/t^ =12.5 for 
aluminium. Similarly as can be seen in figures 3.7 to 3.12, in the 
case of roll torque, except for the case of aluminium with 
R/h.^39, in all cases the deviation between analytical and 
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experimental results Increases. While it is reported that the 
deviation is because of the assumption of constant coefficient of 
friction [ 15 can also be attributed to the fact that at 

large values of reduction it becomes difficult to predict the 
neutral point accurately , which has a direct effect on the roll 
torque. Although the effect may not be at par with that of 
friction, the rigid roll assumption and the ignorance of entry 
and exit elastic effects are also expected to have contributed to 
the deviation of analytical and experimental results. Figures 
3.13 to 3.16 show the comparison of analytical and experimental 
normal pressure distribution along the arc of contact. Although 
there is a considerable variation in magnitude, the Important 
feature of multiple peaks of normal pressure at low values of 
R/h^ and the friction-hill type of pressure distribution at 
higher values have been reproduced. The multiple peaks at low 
R/h^ have been predicted also by the slip-line field solution 
of Flrbank et al [ 8 J. This suggests that the present model 
would be able to predict correct magnitudes of normal pressure 
over some range of values of R/h^ reasonably accurately. 

3.2. Parametric Study 

In this section a discussion on the effect of process 
variables viz. reduction ratio, R/l^ and the coefficient of 
friction on the normal pressure distribution, shear stress 
distribution, roll separating force and the roll torque is 
presented . 

3.2.1. Effect of reduction ratio and R/l^ 

The figure 3.17 depicts the variation of normal pressure(t n ) 
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along the arc of contact at different reductions. While the 
distribution is more uniform at higher reductions, at low 
reductions it is non-uniform and exhibits multiple peaks. It is 
interesting to note the peak point in the normal pressure curve 
at the entry which is the characteristic of all distributions. 

The normal pressure, as shown in figure 3.18, increases 
rapidly with decreasing values of R/h^ and becomes less smooth. 
Again the interesting observation is at low values of R/h 1 there 
is a distinct peak at the entry which is not to be seen at high 
values of R/h 1 , for example, R/h.^130. However, at the exit the 
variation is smooth and uniform for all values of R/t^. 

Figure 3.19 shows the distribution of shear stress along the 

arc of contact from exit to the entry at different reductions. 

Since It |=f|t I, at higher reductions t is more uniform similar 

to t . From figure 3.20, where the distribution of shear stress 
n 

at different values of R/h^ is shown, it is clear that the peak 
values of t g also, similar to t^, increase rapidly with 
decreasing R/h^. 

The variation of roll separating force with R/h^ at 
different reductions Is shown in figure 3.21. The figure suggests 
that while the roll separating force decreases with Increasing 
R/hj^, at all reductions in the range, this decrease becomes more 
rapid at Increasing reductions. This is consistent with the trend 
of t with R/h 1 - This may be attributed to the indenting effect 
of roll into the work piece which increases with decreasing R/l^. 
At a given R/t^, roll separating force increases with reduction 
In spite of decrease in t ft . So the increase in contact area seems 
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to be having more effect than decreasing t on the roll 
separating force. Similar observations can be made in the figure 
3.22 which shows the variation of roll torque with R/l^ at 
different reductions. However, in the case of roll torque the 
decrease with R/h. at higher reductions is more monotonous. The 
continuous increase in roll torque with reduction at any given 
R/h 1 can be justified by referring to the variation of shear 
stress in figure 3.19, where it is clear that the percentage 
length of the backward slip zone (distance from the entry to the 
neutral point ) at any given reduction increases with Increasing 
reduction. Thus it follows that since the shear stress is 
positive in thebackward slip region roll torque increases with 
Increasing reductions. 

3.2.2. Effect of the coefficient of friction 

The figure 3.23 shows the distribution of normal pressure 
along the arc of contact at different coefficients of friction. 
Uhlle the shape of the curve and smoothness are restored for all 
values of coefficient of friction, the peak magnitude increases 
with increasing coefficient of friction and peak seems to shift 
towards the entry. Figure 3.24 which gives the variation of shear 
stress along the arc of contact at different coefficients of 
friction. It can be noted from this figure that the peak 
magnitude increases with increasing coefficient of friction and 
that the neutral point shifts towards entry. 

Figures 3.25 and 3.26 show the variation of roll separating 
force and roll torque with respect to the coefficient of 
friction, respectively. The modes of variation of roll force and 
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roll torque In these graphs can be appreciated by referring to 
the figures 3.23 and 3.24 where the normal pressure and shear 
stress distributions are depicted, respectively. It is 
interesting to see that since the contact area remains the same 
for all the three cases for which the curves are drawn, the 
variations, whatsoever, in roll force and roll torque are brought 
about by the change in shape and mode of variation of normal 
pressure and shear stress curves. From figure 3.25 it is evident 
that the roll force first increases with the coefficient of 
friction, however gradually, and then becomes uniform. This seems 
to be the consequence of the variation of normal pressure 
distributions which seem to become similar bothi^shape and 
magnitude with Increasing coefficients of friction. Further the 
decreasing trend of roll torque curve with increasing coefficient 
of friction after reaching a peak can be appreciated by observing 
the figure 3.24. Here it is very clear to observe that with 
increasing coefficient of friction not only the neutral point 
shifts towards the middle of the arc of contact but the 
difference between negative and positive average values of shear 
stress decreases, both leading to a net decrease in the roll 
torque with increasing coefficient of friction. However, 
meanwhile the fact that the computer program in this analysis 
ceases to converge for the values of coefficient of friction 
beyond 0.29 should also be pointed out. 

3.3. Effect of Strain Hardening 

In this section the effect of strain hardening is brought 
into light by selecting two sets of process variables. In the 
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first set the process variables are such that they are expected 
to give a minimum value of the design parameter to be obser ed 
and in the second set they give an expected maximum of the same. 
These expectations are based on the foregoing parametric study. 
The design parameters are the same which were considered for 
parametric study. Uhile the discussion for the former two is 
given by means of graphs, the later two are presented in the 
table 3.2. The percentage deviation for peak magnitudes of normal 
pressure and shear stress are also included in the table 3.2. 

The figure 3.27 contains the distribution of normal pressure 
without and with strain hardening for the two sets of process 
variables. Uhile it is clear to see that the normal pressure is 
higher with strain hardening than without strain hardening all 
along the arc of contact in both the sets, the peak deviation 
occurs near exit in the first set whereas it occurs in the middle 
of the arc of contact in the second set. Further the percentage 
error in the estimation of t when the strain hardening effect is 
neglected is more for the set 2 (see the table 3.2). It can be 
noted that the variation of along the arc of contact is more 
uniform for the set 1 than for set 2. Figure 3.28 shows the 
effect of strain hardening on the shear stress distribution along 
the arc of contact for the same two sets of process variables. It 
shows that the percentage deviation is more for set 2. The 
underestimation of roll force and roll torque by not considering 
the effect of strain hardening is presented in table 3.2. 
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3.4. The Study of Contours 

The figures 3.29 to 3.32 show the contours for the 
normal-strain rate, equivalent strain-rate, deviatorlc stress 
invariant and the equivalent strain in the domain for four sets 
of values of the three process variables. The four sets are so 
selected that the comparison with the contours of the first set 
of any of the successive three sets gives the effect of one 
process variable on the distribution of that parameter in the 
domain. The graphics routine for plotting these contours has been 
developed along with the computer program for the main finite 
element analysis. 

The figure 3.29 shows the contours for normal strain rate 
for steel for the four sets of process variables in (a),(b),(c) 
and (d). The normal strain rates, at any given set of process 
variables first increase from a small value as the metal enters 
the deformation zone and reach a maximum and thereafter decrease 
gradually as the metal exits from the deformation zone. However, 
a marked difference can be observed between these trends along 
the surface where the work piece is in contact with the roll and 
along the axis of symmetry. While the normal strain rates along 
the axis of symmetry uniformly increase, reach a maximum and then 
decrease to a small value, just as explained above as a typical 
nature, along the work piece-roll interface they undergo a 
non-uniform variation with multiple maxlmums and minimums from 
inlet to the exit. This variation from uniform along the axis of 
symmetry to non-uniform along the interface is gradual. Because 
of this inhomogeinity of deformation pockets of normal strain 
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r At£8 rfiflult Along and just below the lntfirfscc . It c&n be 
observed from the figure 3.29 (a), (b), (c) and (d) that with 
increased reduction ratio ( a to b ) the homogeinity of 
deformation increases and with increased R/h^ and coefficient of 
friction the inhomogeinity increases ( a to c and a to d, 
respectively ). The normal strain rates in 3.29(a) first 
increase and reach a maximum value at the first contact point. 
Then it decreases towards the middle of the arc of contact 
reaching a minimum. It then again increases and attains another 
peak value and continuously decreases towards exit. They have two 
peaks on the surface and a minimum, Uith increased reduction, as 
shown in 3.29(b), the number of peaks reduce to one which occurs 
at a point away from the first contact point. Further the normal 
strain rates continuously increase till the peak is reached and 
then continuously decrease till they reach the cut off value at 
the exit. The Increasing coefficient of friction and decreasing 
R/h^ seem to have little effect on the distribution of normal 
strain rates, as can be seen from 3.29(c) and 3.29(d). However, a 
minimum Is Introduced before a contour goes down to cut off value 
at the exit. Similar observations are possible from the 
equivalent strain rate contours shown in the figure 3.30. 

The contours of deviatorlc stress Invariant which is a 
measure of stress at any point required to cause a given amount 
of plastic strain at that point . In general, at any given set of 
process variables, the equivalent deviatorlc stress first rapidly 
increases at the entry as the loading occurs, continues to 
increase Into the deformation zone reaching a maximum somewhere 
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In the middle and then gradually decreases towards exit. However, 
at the exit where unloading of the work piece occurs , the 
decrease In the devlatoric stress invariant Is again very rapid 
as the metal comes out of the deformation zone. Although the 
effective devlatoric stress remains more or less uniform across 
the rolling direction in the work piece at the loading and 
unloading points the variation between and surface particles is 
distinct. This observation is common to all sets in the figure 
3.31 except that the variations at the inlet and exit are not 

shown in any of these figures. Comparison of figure 3.31(b) with 
figure 3.31(a) suggests that the devlatoric stresses become more 
or less uniform across the section with Increased reductions. 
However, there is not much effect of other two process variables 
on the distribution of these contours. The figure 3.32 shows the 
contours of equivalent strain, also in four sets. The following 
observations can be made from the figures for the contours of 
equivalent strain. At the entry as the metal enters into to the 

deformation zone the strain increases rapidly from a very low 

value. In this zone clearly the surface particles undergo higher 
amounts of strain than those Inside. As the metal progresses 
towards the exit, however, the Inside particles catch up and 

trend becomes almost reverse. While the trend at higher reduction 
ratios is more or less similar, because of longer path available 
for the particles to travel from the inlet to the exit this 
variation Is gradual. Equivalent strains at higher coefficient of 
friction are shown in figure 3.32(a) keeping the other variables 
constant. This illustrates that at higher friction the surface 
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partlclsa &re too highly strained to be reached by the inside 
particles and therefore throughout the domain the equivalent 
strain across the sections increases from inside to the surface. 
This variation is gradual from a low value to a high value. 

3.5. Conclusions 

Based on the above discussions the following conclusions can 
be drawn. In the present analysis cold plane-strain flat rolling 
problem is solved by the finite element method at steady state by 
considering the material to be rigid-plastic strain hardening. 
The model can give good results over a range of the process 
variables viz. the reduction ratio, R/h^ and the average 
coefficient of friction. The model has the capability of 
accommodating variable coefficient of friction which can be 
considered from the experimental results to obtain better 
results. In this work a detailed parametric study has been 
conducted to manifest the effect of the process variables on the 
Important design parameters of rolling process. The error of the 
assumption of perfectly plastic material is highlighted. Contours 
have been shown to illustrate the distribution of stresses, 
strain rates and strains for an understanding of the physical 
process of rolling. 

3.7. Scope for Further Research 

As has been indicated in the above lines, the model is 
certainly amenable to the inclusion of the variable coefficient 
of friction which can further improve the accuracy of results. It 
should also be pointed out that, a change in the mesh by bringing 
in smaller elements at appropriate places is likely to improve 
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the results. Uhl 1 e three dimensional treatment Is not an 
improvement over every 2-D analysls( especially when the 
width/ thl ckness is very high the cost incurred in 3-D analysis 
can prove to be too high to appreciate the extra accuracy 
achieved), by carrying 3-D analysis some additional features like 
profile of the strip can be brought into light and the problems 
of shape rolling where the flow of metal in the third direction 
is significant to achieve the required shapes can be dealt with. 
The assumption that the metal is rigid at the entry and exit and 
that the roll is rigid may also be eliminated to predict better 
results . 
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Metal; Steel — Ana. 



Fig. 3-1 Comparison ot Analytical and Experimental 
Roll Force 


Metal' Steel — Ana. 



Fig. 3-2 Comparison ot Analytical and Experimental 
Roll Force 
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Metal: Aluminium — Ana. 

t =0‘10,R/h-| =12*5 .h^^Amm □ Exp. [25] 


Reduction (%>) 

Fig. 3-3 Comparison ot Analytical and Experimental 
Roll Force 


Metal : Aluminium — Ana. 

f =0*10, R /h-) =39,^=2*032 mm □ Exp. [2.5] 


? 5 


■o 

$ 3 



O 1 


0 10 30 50 

Reduction (°/o) 

Fig. 3-4 Comparison ot Analytical and Experimental 
Roll Force 






Fig. 3-6 Comparison of Analytical and Experimental 
Roll Force 
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Fig.3-8 Comparison ot Analytical and Experimental 
Roll Torque 
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Metal : Aluminium Ana. 

E t = 0-10, R/h 1 =125, ^=6-276 mm a Exp.[25j 



Fig.3-9 Comparison of Analytical and Experimental 
Roll Torque 


Metal : Aluminium — Ana. 



Fig.3-10 Comparison ot Analytical and Experimental 
Roll Torque 
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_ Metal -.Copper — Ana. 

£ t=0-06,R/h-p12-5,h-|=6-35mm □ Exp. [25] 



Fig.3-11 Comparison of Analytical and Experimental 
Roll Torque 



Fig. 3-12 Comparison of Analytical and Experimental 
Roll Torque 
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10 Normal Pressure (GN/sq m) 

Normal Pressure (GN/sq.m) w ^ 


Metal : Aluminium ---Ana. 

t = 0-10, R/h-j =12*5, b|=6-274,7 0 r =14*17— Exp. [25J 



•13 Comparison of Analytical and Experimental 
Normal Pressure 


Metal’- Aluminium — Ana. 

t = 0-10, R/1^39, 1^=2032,7. r s 33-75 —Exp. [25] 



Fig. 3-14 Comparison ot Analytical and Experimental 
Normal Pressure 
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Arc of contact (degrees) 


Fig.3-15 Comparison of Analytical and Experimental 
Normal Pressure 


MetaF Copper --Ana. 

i=0-06,R/h 1 =50,h 1 =1-562,%r=22-76 — Exp.[2S] 



Fig.3-16 Comparison ot Analytical and Experimental 
Normal Pressure 
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Metal : Steel 8 °/ c 

f = 0.06, R/h|= 65, R = 65 mm 16°/ 
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Fig.3J7Distribution of normal pressure at different reductions 
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(Exit) (Entry) 

Fig. 3.19 Distribution of shear stress at different reductions. 
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R/h| 

Fig. 3.21 Variation of roll separating force with R/h t at d 
reductions. 



Roll torque/Unit width ( kN-m/m) 


Metal : Steel f = 0.06, R = 65 mm. 



Fig. 3. 22 Comparision of variation of roll torque with R/h( at 
different reductions. 
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Fig.3.23 Distribution of normal pressure at different coefficients of friction 
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Fig3.27 Distribution of normal pressure with and without strain hardening 
(a) minimum and (b) maximum in the range'. 
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Fig. 3.28 Distribution of shear stress with and without strain hardening (a) minimum and 
(b) maximum in the range 





Tabl* 1 . Material Properties and Stress-Strain Relations 


Material 

Density 

In the Relation, <y y =o- o (l+K2) n 


■MM 


Ck*/* 3 ) 

<y Q C MPa ) 

n 

1/ K 

Reference 


Aluminium 

2800 

50.3 

0.260 

0.050 

Al-Salehi [ 25 

] 

St eel 

7870 

358.0 

0.300 

0.044 

Shida [ 26 ] 


Steel 

7870 

324.0 

0.295 

0.052 

Shida [ 26 ] 


Copper 

8600 

70.3 

0.490 

0.022 

Al-Salehi [ 25 

] 


Table 2. The Percentage underestimation in the analysis without 
strain hardening for roll separating force, roll torque, peak 
magnitude of normal pressure and absolute peak magnitude of shear 
stress 


Parameter 

Set 


Parameter Value 

% Error 

(w.r.t. the value 
with strain hard.) 




— . 

without 

with 

Roll Force 
CkN/m) 

Set 

1 

596 

707 

15.70 

Set 

2 

2075 

3111 

33.45 

Roll Torque 
(kN-m/m) 

Set 

1 

1.06 

1.15 

07.82 

Set 

2 

11.13 

16.35 

31.93 

(t ) 
n max 

( MPa ) 

Set 

1 

382 

485 

21.24 

Set 

2 

700- 

1160 

39.66 

I<1 

( MPa ) 

Set 

1 

24 

30 

20.00 

Set 

2 

135 

230 

41.30 
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* no. of nodes 261 

* no. of elaents 56 




(b) *r = 32 R/h 1 = 130 R = 65 mm f = 0.06 



(c) *r = 8 R/hj =65 R = 65 am f = 0.06 



Fl fi . 3.2? Normal Strain-Rate Contoure at Four Different Seta of 
Proceaa Variablea, (a), (b), (c) and (d). 
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* no. of nodes 261 

* na of elements 56 
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